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Vector decay constants in quarkonia 
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Lattice NRQCD with leading finite lattice spacing errors removed is used to simulate heavy-heavy vector decay 
constants. Quenched simulations are performed at three values of the coupling and fifteen values of the quark mass. 
The improved gauge action with plaquettes and rectangles is used. Landau link improvement is used throughout. 
"Perturbative" and nonperturbative meson masses are compared. One-loop perturbative matching between lattice 
and continuum heavy-heavy vector currents is performed. The data is consistent with afv oc yJMv a. 



1. INTRODUCTION 

The gross features of heavy quarkonia are well 
described by lattice NRQCD [1-3]. However, spin 
splittings (which are small for quarkonia) tend 
to be underestimated — even when relativistic cor- 
rections are included Q|. Spin splittings are one 
measure of the mesonic wave function at the ori- 
gin; another is the vector meson decay constant. 

A previous study || °f t nc vector decay con- 
stant found large corrections from the perturba- 
tive matching between the continuum and lattice 
matrix elements. Another study reported a 
rather imprecise value for simulations with the 
order v 2 classically improved vector current. In 
this paper, we try to improve these calculations 
by removing the leading finite lattice spacing er- 
rors in the fermion and gauge actions in a sym- 
metric fashion and by performing a more precise 
simulation with the inclusion of the order v 2 clas- 
sically improved current. 

1.1. Lattice NRQCD 

The fermion Lagrangian is discretized in a sym- 
metric fashion with the leading finite lattice spac- 
ing errors in the spatial and temporal derivatives 
removed — all links are tadpole improved by divid- 
ing by uqi the average link in the Landau gauge: 
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n is the stability parameter chosen to satisfy 
n > — . A^ 2 ) is the gauge-covariant lattice 

ma G G 

Laplacian, and A^ 4 ) is the gauge-covariant lattice 
quartic operator (^\ Df). 

The gauge action is tadpole improved with 
leading finite lattice spacing errors removed by 
rectangles: 
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The 'kinetic' meson mass is defined by 



E(p) - E(0) = 
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where E{p) is the simulated meson energy. 

1.2. Perturbative lattice NRQCD 

The Feynman rules are derived from the La- 
grangians of Eqs. ([!]) and (|3|) by making the re- 
placement U/j,(x) — * exp[iagA'^(x)T a ]. The gluon 
propagator follows from the quadratic piece of the 
gauge action in Eq. (||) — we take the piece pro- 
portional to <5 p v in this paper. 

The continuum and lattice decay constants 
in one-loop perturbation theory are related by 

fy = Z match fv,latt, where Z ma tch = 1 + 

g 2 (— gfy - {SV + SZ)i att ). Note that a linear in- 
frared divergence cancels between the continuum 
and lattice perturbative shifts. Fig. [j] shows the 
results of the one-loop calculation. The shown er- 
rors are estimates of the extrapolation errors. All 
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Table 1 

Mp er t versus Mkm for the runs nearest the charm and bottom regions respectively. 

Mp er ta 
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a(fm) 


ma[n] 


MkinCL 


w/o tad imp 


w tad imp 


M fem (GeV) 


7.2 


.240(3) 


1.6[3] 


3.37(6) 


2.89(7) 


3.60(7) 


2.77(6) 
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1.5[3] 


3.20(5) 
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3.37(6) 


3.08(6) 
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.178(3) 


1.4[3] 


3.00(6) 
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3.15(6) 
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.198(2) 


4.4[2] 


9.07(21) 


9.47(68) 


8.98(68) 


9.05(23) 


7.3 


.174(3) 


4.3[2] 


8.74(20) 


9.17(61) 


8.74(61) 


9.91(28) 
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3.5[2] 
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Figure 1. Perturbative lattice NRQCD one-loop 
contribution to the matching of the heavy-heavy 
vector current. These results must be multiplied 
by g 2 - 

one-loop integrals are performed by VEGAS 
in this paper. 

The "perturbative" meson mass is defined by 
M pe rt = 2(mZ m - -B ) + E slm . Table | shows 
the simulation parameters near the physical re- 
gions and compares M per t and Mkin- The lattice 
spacing is fixed by the S-P splitting: 458 MeV, 
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and 



a boosted coupling is used: g 
uq = 1 — j^u^ is used to define u^. 

1.3. Decay constants 

The asymptotic form of a meson propagator is 
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where J(x,t) = XxF x ip x is a non-relativistic 
current with T x = Q x ^/ X , where VL X interpo- 



lates the meson of interest and 7^ is a smear- 
ing operator chosen in a gauge-invariant fashion: 
7^ = [l + £A< 2 )(i)]"' . We set e = 1/12 and tune 
the smearing parameter n s to maximize the over- 
lap with the state of interest. We find the range 
7-30 for n s to be sufficient, with the P-wave re- 
quiring about twice as much smearing as the S- 
wave and the smearing parameter increasing for 
decreasing quark mass. 

The above matrix element for a vector at 
rest is related to its decay constant fy through 
(0|J(0)|F) = ^7=, where a non-relativistic 
norm is used. 

Finally, note that simulations are performed 
with an order v 2 classically improved current. 
The improved interpolating operator for a non- 
relativistic vector meson is given by 
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where A is the symmetric lattice gauge-covariant 
derivative. 

2. RESULTS AND DISCUSSION 

The data sample includes 1600 quenched gauge 
field configurations at /? = 7.2 and 7.3 (8 3 x 14), 
and 1200 configurations at = 7.4 (10 3 x 16). 
Multiple sources along the spatial diagonal are 
used to measure the local-smeared meson corre- 
lators with the number of smearing steps opti- 
mized to have the best overlap with the state of 
interest. All plots include bootstrap errors using 
twice as many bootstrap ensembles as there are 
configurations. 



3 




. 15 



. 17 



.19 .2 1 

a ( fmn ) 



.23 



.25 



Figure 2. Scaling behavior of the charm decay 
constants including the perturbative matching. 
The dotted line is the experimental result. 
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Figure 3. Scaling behavior of the bottom decay 
constants including the perturbative matching. 
The dotted line is the experimental result. 



Figs. [2|-[l| show our main results. Note the 10% 
underestimations (order v 2 improved data) as ex- 
pected from the previous work on the spin split- 
tings. Also note the approximate \JMy a depen- 
dence. Shortly after the discovery of charm, Yen- 
nie |^] noticed this same dependence from empir- 
ical data: 
T v - 

-f- ~ constant - 12 keV (7) 

for light through heavy ground-state vector 
mesons (e q is the quark charge in units of e). 
This is the same as our data since the leptonic 
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Figure 4. Meson mass dependence of the order v 2 
improved decay constant including the perturba- 
tive matching. The dotted line is proportional to 

y/Mya. 

width is proportional to e 2 j^. A final note is 
that a linear (Coulomb) potential implies a con- 
stant (linear) dependence on the meson mass for 
the decay constant, so our data is consistent with 
a superposition of the two potentials. 
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